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The origin of cosmic dust is a fundamental issue in planetary science. This paper revisits the origin of dust 
in galaxies, in particular, in the Milky Way, by using a chemical evolution model of a galaxy composed of stars, 
intestellar medium, metals (elements heavier than helium), and dust. We start from a review of time-evolutionary 
equations of the four components, and then, we present simple recipes for the stellar remnant mass and yields 
of metal and dust based on models of stellar nucleosynthesis and dust formation. After calibrating some model 
parameters with the data from the solar neighbourhood, we have confirmed a shortage of the stellar dust production 
rate relative to the dust destruction rate by supernovae if the destruction efficiency suggested by theoretical works 
is correct. If the dust mass growth by material accretion in molecular clouds is active, the observed dust amount in 
the solar neighbourhood is reproduced. We present a clear analytic explanation of the mechanism for determining 
dust content in galaxies after the activation of accretion growth: a balance between accretion growth and supernova 
destruction. Thus, the dust content is independent of the uncertainty of the stellar dust yield after the growth 
activation. The timing of the activation is determined by a critical metal mass fraction which depends on the growth 
and destruction efficiencies. The solar system formation seems to have occured well after the activation and plenty 
of dust would have existed in the proto-solar nebula. 
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1. Introduction 

Cosmic dust grains are negligible in mass in the Universe. 
Nevertheless, they play significant roles on a lot of astro- 
nomical, astrophysical, and astrochemical aspects: extinc- 
tion (absorption and scattering) matter of radiation, an emis- 
sion source in infrared wavelengths, a coolant and a heat 
source in the interstellar medium (ISM) and intergalactic 
medium (IGM), and a site of formation of molecules. There- 
fore, dust is one of the most important ingredients in the Uni- 
verse. Dust is also important for planetary science because 
grains are material for planets. 

Dust grains are formed in rapidly cooling gas of stellar 
outflows (Yamamoto & Hasegawa 1977; Draine & Salpeter 
1977). We call such grains 'Stardust'. Sources of the Stardust 
are asymptotic giant branch (AGB) stars, supernovae (SNe), 
red supergiants, novae, Wolf-Rayet stars, and so on (e.g., 
Gehrz 1989). The main source of the Stardust in the present 
Milky Way and the Magellanic Clouds is thought to be AGB 
stars (Gehrz 1989; Draine 2009; Matsuura et al. 2009). 

SNe may also produce a significant amount of Stardust 
(Kozasa & Hasegawa 1987; Todini & Ferrara 2001; Nozawa 
et al. 2003, 2007; Schneider et al. 2004; see also Kozasa 
et al. 2009). Stardust from SNe is particularly important in 
the early Universe because the time for stars to evolve to the 
AGB phase is typically about 1 Gyr, but the cosmic time in 
the early Universe is shorter than it (Morgan & Edmunds 
2003; Maiolino et al. 2004; Dwek et al. 2007; but see also 
Valiante et al. 2009). The 'first' Stardust may also play an 
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important role to change the mode of star formation from 
massive star dominated to present-day Sun-like star domi- 
nated (Schneider et al. 2003, 2006). 

However, dust formation by SNe remains in controversy 
observationally. First detections of a few M© dust freshly 
formed, which is much larger than expected, in Cassiopeia 
A (Cas A) and Kepler SN remnants (SNRs) by submilimeter 
observations with SCUBA (Dunne et al. 2003; Morgan et 
al. 2003) were almost contaminated by foreground dust in 
the ISM on the sight-lines (Krouze et al. 2004; Gomez et 
al. 2009). Recent infrared observations with Spiter Space 
Telescope and AKARI and submilimeter observations with 
Herschel and BLAST of Cas A and other SNRs are agreeing 
with theoretical expectations of 0.01-0.1 Mq per one SN 
(Rho et al. 2008; Sakon et al. 2009; Nozawa et al. 2010; 
Barlow et al. 2010; Sibthorpe et al. 2010). 

Once Stardust grains are injected into the ISM, they are 
processed there. The grains in hot gas are bombarded by 
thermally moving protons and sputtered (Onaka & Kamijo 
1978; Draine & Salpeter 1979). SN shock waves proba- 
bly destroy dust grains by grain-grain collisional shattering 
as well as sputtering (e.g., Dwek & Arendt 1992; Jones et 
al. 1994, 1996; Nozawa et al. 2006; Silvia et al. 2010). This 
destruction process is widely accepted and observational ev- 
idences of the destruction have been found in several SNRs, 
especially with Spitzer Space Telescope (Arendt et al. 1991, 
2010; Borkowski et al. 2006; WilUams et al. 2006; Dwek 
et al. 2008; Sankrit et al. 2010; but see Mouri & Taniguchi 
2000). 

Assuming the destruction efficiency expected by theoreti- 
cal works, we obtain the life-time of dust grains of the order 
of 100 Myr (McKee 1989; Draine 1990; Jones et al. 1994, 
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1996). On the other hand, the injection time of Stardust is of 
the order of 1 Gyr (e.g., Gehrz 1989). Thus, another efficient 
channel of dust formation is required to keep dust content 
in galaxies. The most plausible mechanism is the accretion 
growth in the ISM (Draine 1990, 2009); in dense molec- 
ular clouds, atoms and molecules of some refractory ele- 
ments and compounds accrete onto pre-existent grains and 
may change from the gas phase to the solid phase. Note that 
unlike the sticking growth of grains well studied in proto- 
planetary disks, this accretion growth causes an increase in 
dust mass. This type of growth is favored to explain the ob- 
served depletions of some elements in the gas phase of the 
ISM relative to the solar abundance. The correlation between 
the depletion degree and the density in the ISM is particularly 
suggestive for this process (e.g.. Savage & Sembach 1996; 
Jenkins 2009). It is also suggested that an efficient growth is 
required to explain massive dust mass observed in the early 
Universe (Michalowski et al. 2010). 

Since the pioneering work by Dwek & Scalo (1980), 
many theoretical works on dust content evolution in galax- 
ies have been made so far (Dwek 1998; Lisenfeld & Fer- 
rara 1998; Edmunds & Bales 1998; Hirashita 1999a,b,c; 
Edmunds 2001; Hirashita et al. 2002; Inoue 2003; Mor- 
gan & Edmunds 2003; Dwek et al. 2007; Zhukovska et 
al. 2008; Calura et al. 2008; VaHante et al. 2009; Pipino 
et al. 2011; Dwek & Cherchneff 2011; Gall et al. 2011a,b; 
Mattsson 2011; Asano et al. 2011). These works are based 
on the evolutionary model of elemental abundance in galax- 
ies called chemical evolution model (Tinsley 1980 for a re- 
view) and add some (or all) of the three processes of forma- 
tion, destruction, and growth of dust to it. One of the main 
results from the recent works is the importance of the accre- 
tion growth. 

This paper presents a new interpretation of the mechanism 
for determining dust content in galaxies. Previous works im- 
ply that the mechanism is a balance between dust destruction 
by SNe and accretion growth in the ISM. However, this point 
has not been discussed clearly, in contrast, this paper analyt- 
ically shows that it is. For this aim, a simple one-zone model 
is sufficient. In addition, we present new simple recipes de- 
scribing stellar remnant mass and yields of elements and dust 
from state-of-the-art models of stellar nucleosynthesis and 
formation of Stardust. 

We will start from a review of basic equations presented in 
§2. In §3, we present new simple recipes of stellar remnant 
mass and yields. In §4, we calibrate some model parameters 
to reproduce the observed properties of the solar neighbor- 
hood. In §5, we present our analytical interpretation of the 
mechanism for determining dust content in galaxies. We will 
present some further discussions in §6. Experts of this field 
may go straight to §5 which is the new result of this paper. 

Throughout this paper, we call elements heavier than he- 
lium 'metal' according to the custom of astronomy. We 
adopt the metal mass fraction (so-called metallicity) in the 
Sun of Zq — 0.02 (Anders & Grevesse 1989) convention- 
ally, although the recent measurements suggest a smaller 
value of 0.0134 (Asplund et al. 2009). 



2. Chemical and dust evolution model of galaxies 
2.1 Equations of chemical and dust amount evolution 

We deal with a galaxy composed of stars (including their 
remnants; i.e. white dwarfs, neutron stars, and black-holes) 
and the ISM. For simplicity, we assume the ISM to be one- 
zone. The ISM contains metal and dust as internal com- 
ponents. If we denote masses of these components as Af, 
(stars [and remnants]), Mjsm (ISM), Mz (metal), and 
(dust), the equations describing their time evolutions are 
(e.g., Dwek 1998) 



dt 



= S{t) - R{t) , 



dNh 



ISM 



dt 



dMz 
dt 



= -S{t) + R{t) + I{t) ~ 0{t) , 
Z{t)S{t) + Yz{t) + Iz{t) - Oz{t) . 



(1) 



(2) 



(3) 



dt 



= -Z^{t)S{t)+Yd{t)^Ds^{t)+G, 



.{t)+Idit)-Od{t) , 
(4) 



where S is the star formation rate, R is the mass return rate 
from dying stars, Yz and Id are the metal and dust supplying 
rates 'yields' by dying stars, respectively. Z = Mz/Mism 
is the metal mass fraction in the ISM called 'metallicity', and 
Zd = Md/MisM is the dust mass fraction in the ISM which 
we call the dust-to-gas mass ratio. Note that A/ism > AIz > 
A/d. 

/, Iz, and Id are the ISM, metal, and dust infall rates from 
the IGM, respectively. O, Oz, and Od are the ISM, metal, 
and dust outflow rates to the IGM, respectively. In this paper, 
we do not consider any outflows (O ~ Oz = Od = 0), 
but consider only an ISM infall / (no metal and dust in 
infalling gas: Iz = Id = 0), which is required to reproduce 
the metallicity distribution of stars nearby the SunQ The 
reason why we omit any outflows is that we do not know the 
transport mechanism of metal and dust from galaxies to the 
IGM (e.g., Bianchi & Ferrara 2005). However, this omission 
may be inconsistent with detections of metal and dust in the 
IGM (e.g., Songaila & Cowie 1996; Menard et al. 2010)0 

In the dust mass equation (eq. [4]), there are two additional 
terms; Dsn is the dust destruction rate by SNe and Gac is the 
dust growth rate in the ISM by metal accretion. These two 
terms are discussed in §2.5 and §2.6 in detail. 
2.2 Star formation and infall rates 

We adopt a simple recipe for star formation introduced by 
Schmidt (1959): S oc Mism^ (Schmidt law). The index 
p called Schmidt index is observationally indicated to be 



^Without gas infall from the intergalactic space, we expect a much larger 
number of low-metallicity stars at the solar neighborhood than the observed. 
This is called the 'G-dwarf problem' (e.g., Pagel 1989). 

^The origin of intergalactic metals and dust is galactic outflows and the 
amount ejected from galaxies is the same order of that remained in galaxies 
(e.g., Menard et al. 2010 for dust; see also Inoue & Kamaya 2003, 2004, 
and 2010). Dust grains may be ejected from galaxies more efficiently than 
metals because the grains receive momentum through radiation pressure 
(Bianchi & Ferrara 2005). Even in this case, our discussion about the dust- 
to-metal ratio in §5 would not be affected essentially by omission of this 
selective removal of dust, although the set of model parameters which can 
reproduce the observations would change. In any case, this point would be 
an interesting future work. 
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p = 1-2 (e.g., Kennicutt 1998; Elmegreen 2011) and some 
theoretical interpretations for the value are presented (e.g., 
Dopita & Ryder 1994). However, the value and its origin 
of the index is still an open problem (Elmegreen 2011 and 
references therein). Fortunately, the choice of the index is 
not important in fact because in §4 we calibrate other model 
parameters so as to reproduce the observed star formation 
history S{t) at the solar neighborhood which is essential. 
We here assume p — 1 in order to solve the equations 
analytically in §5. In this case, we need a time-scale to give 
the star formation rate: star formation time-scale, tsf (see 
Table 1 in §4 for the values). Thus, the star formation rate is 
given by 



S{t) 



Tsf 



(5) 



The infall from the IGM mimics the structure formation 
in the Universe based on the hierarchical scenario with cold 
dark matter (e.g.. Peacock 1999); small galaxies are first 
formed at density peaks of the dark matter distribution in the 
Universe and they grow up larger and larger as they merge 
each other and also obtain mass by accretion process. Here 
we simply assume a smooth exponential infall rate although 
the mass assembly of a galaxy is intrinsically episodic due 
to the merging process. This simplification is a kind of en- 
semble average of many galaxies and appropriate to examine 
a mean property of the galaxies. The infall rate which we 
adopt is 

Mtotal 



m 



■exp(-t/rin) . 



(6) 



where Tin is the infall time-scale and Mtotai is the total mass 
which a galaxy obtains within the infinite time (see Table 
1 in §4 for the values). Note that Mtotai just gives the 
normalization of mass of a galaxy. 
2.3 Stellar mass spectrum and returned mass rate 

Salpeter (1955) first investigated the mass spectrum of 
stars in the solar neighborhood, corrected it for the modu- 
lation by stellar evolution and death, and obtained the mass 
spectrum of stars when they are bom, called initial mass 
function (IMF) of stars. The Salpeter's IMF is a power-law: 
dN/dm = (/)(m) cx m""^ with q = 2.35. A lot of following 
researches confirmed that the slope was quite universal, es- 
pecially for massive stars, although there was a cut-off mass 
for low mass stars (e.g., Kroupa 2002, Chabrier 2003 for re- 
views). We adopt here a simple functional from proposed by 
Larson (1998) which is essentially equivalent to the IMFs by 
Kroupa (2002) and Chabrier (2003) as 

(f){m) oc exp(— mc/m) , (7) 

with a cut-off mass rric and the range from miow to mup. 
As a standard case, we adopt p = 2.35, nic = 0.2 Mq, 
m\ovj =0.1 Mq, and ruup = 100 Mq. The cut-off mass 
well matches with the observed data compiled by Kroupa 
(2002). We normalize the IMF as J^^^"'' m(f>{m)d'm = 1. 
The mass returned rate from dying stars, R, is given by 

R{t)= {m-w{m,Z[t'])}(j){m)S{t')dm, (8) 

Jm\f(t) 



where 



is the time at which stars with mass m dying at time t are 
bom, Tif (m) is the stellar life-time, w{m, Z) is the remnant 
mass of stars with mass m and metallicity Z, and m\i{t) 
is the minimum mass of stars dying at time t. This is the 
inverse function oft = T\i{m). If time t is less than the life- 
time of the star with mup, the returned rate i? = 0. We have 
assumed that the metallicity of a star is the same as the ISM 
metallicity at the time when the star is bom. 

The stellar life-time Tit(m) is calculated by the formula 
of Raiteri et al. (1996) which is a fitting function of Padova 
stellar evolutionary tracks (Bertelli et al. 1994). This formula 
is a function of stellar mass m and metallicity Z. However, 
the Z-dependence is weak. Thus, we neglect it (we always 
set Z = Zq in the formula). 
2.4 Stellar yields of 'metal' and dust 

When stars die, they eject substantial mass of metal and 
dust into the ISM. The term driving the time evolution of 
metal mass given by equation (3) is the metal supplying rate, 
Yz, called metal yield. Using the IMF, 0(m), and the star 
formation rate, S{t), we can express the metal yield as 



J 77 



mz{m,Z[t'])(l){m)S{t')dm, (10) 



mit(t) 

where m,z is the metal mass ejected from a star with mass m 
and metallicity Z, and t' is given by equation (9). 

The dust supplying rate. Id. called dust yield can be ex- 
pressed likewise: 

Yd{t)= md{m,Z[t'])^{m)S{t')dm, (11) 

where md is the dust mass ejected from a star with mass m 
and metallicity Z, and t' is given by equation (9). 
2.5 Dust destruction by supernova blast waves 

Dust grains are destroyed by SN shock waves due to shat- 
tering and sputtering (e.g., Dwek & Arendt 1992). This dust 
destruction are observed in some SNRs as described in § 1 . 
In this paper, we adopt the dust destruction rate by SNe pro- 
posed by Dwek & Scalo (1980) and McKee (1989): 



TSN(i) ' 



(12) 



where the destraction time-scale tsn is defined as the time- 
scale which all the ISM is swept by 'dust destractive' shock 
waves: 

MismW 



tsn(0 = 



(13) 



emsN-RsN(i) ' 

where i?sN is the SN occurrence rate, msN is the mass swept 
by a single SN, and e is the efficiency of the dust destraction. 
The SN occurrence rate is given by 



/.4OM0 

RsN{t) = / <l){m)S{t')dm , 



(14) 



t' = t — T\f{m) 



(9) 



where we have assumed the mass range for SNe to be 8- 
40 Mq (Heager et al. 2003) and f/ is in equation (9). If 
t < Tif (40 Mq), i?sN = 0. Note that we consider only Type 
II SNe and neglect Type la SNe. The reason is discussed in 
§3. 
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The effective mass swept by dust destructive shock wave, 
emsN is the important parameter. It is estimated to be 
- 1000 Mq, namely e - 0.1 and tosn ~ 10"* M© (Mc- 
Kee 1989, Nozawa et al. 2006). Recent models for starburst 
galaxies in the early Universe often assume an effective mass 
of emsN ^100 Mq which is a factor of 10 smaller than our 
fiducial value (Dwek et al. 2007; Pipino et al. 2011; Gall 
et al. 2011a). Their argument is that starburst activity pro- 
duces multiple SNe which make the ISM highly inhomoge- 
neous and the dust destruction efficiency decreases in such 
medium. However, the solar neighborhood is not the case, 
and thus, we keep ernsN ^ 1000 Mq. 
2.6 Dust growth by 'metal' accretion in the ISM 

In the ISM, atoms of some refractory elements (or refrac- 
tory molecules) may accrete onto a dust grain and may be- 
come a part of the grain. We call this process the accretion 
growth of dust in the ISM (Draine 1990). Note that this pro- 
cess does not need nucleation, and thus, can occur even in the 
ISM. A simple estimate of the growth rate is (e.g., Hirashita 
2000) 



Gac(t) = X,oldNd{t)7Ta^SzVzpf'it) , 



(15) 



where Xcow-^d is the number of dust grains in cold dense 
clouds, a is the grain radius, sz is the sticking probability 
of accreting metals (atoms or molecules), vz is the thermal 
velocity of the accreting metals and p^^^ is the mass density 
of the accreting metals in the gas-phase. Note that all the 
quantities except for X^oidNd in equation (15) are typical (or 
effective) values averaged over various grain radii, elements, 
and ISM phases. The gas-phase metal density is reduced to 
= p\fy^Z{l — 5), where pfgM is an effective ISM mass 
density. We define it as a mass-weighted average density 
of various ISM phases, and then, it is determined by the 
density of dense molecular clouds where the dust growth 
occurs. Note that 5 = M^/Mz, the dust-to-metal mass ratio 
(the dust depletion factor is 1 — (5). For spherical grains, 
Nd = 3Md/(47ra^c7), where a is the typical material density 
of grains. 

Equation (15) can be reduced to 



Gac(t) = 



Tac(t) 



The accretion growth time-scale t^c is 



T'acO 



Z{t){l-S[t]) 



(16) 



(17) 



where the normaUzation rac,o is the parameter determining 
the process: 



T'ac.O 



4ac7 



SXcoldSzVzPfsM 



(18) 



This time-scale is very uncertain, but we will obtain 
Tacfi = 3 X 10^ yr as the fiducial value in §4.2 in order 
to reproduce the dust-to-metal ratio at the solar neighbor- 
hood with the SN destruction efficiency of emsN = 1000 
Mq. This value can be obtained with a set of parameters 
of a = 0.1 iJ,m (typical size in the ISM of the Milky Way), 
C7 = 3 g cm~^ (compact siUcates), sz = 1, = 0.2 km s~^ 



(°°Fe as an accreting metal atom and thermal temperature of 
100 K), pfl^ = 1 X 10-22 g cm-3, and X^^id = 0.2. This 
set is just an example but ensures that the time-scale is not 
outrageous. 

There is a discussion about the lifetime of dense clouds 
(or recycling time-scale of dense gas) should be longer than 
the accretion growth time-scale for an efficient dust growth 
(Zhukovska et al. 2008; Dwek & Cherchneff 201 1). Accord- 
ing to these authors, the lifetime is long enough to realize an 
efficient dust growth in the Milky Way and even in starburst 
in the early Universe. Another issue is the effect of grain size 
distribution which is discussed in Hirashita (201 1). 

3. Stellar remnant and 'metal' and dust yields 

In this section, we present new simple formulas to describe 
the stellar remnant mass and yields of metal and dust which 
are useful to input into chemical evolution codes. We rep- 
resent all elements heavier than helium as just a 'metal' in 
the formulas for simplicity, while yields of various elements 
are presented in literature. We consider three types of stellar 
death: white dwarfs through the AGB phase, core-collapse 
Type II SNe, and direct collapse toward black-hole called 
'collapser' (Heger et al. 2003). In this paper, we assume the 
mass range for the SNe to be 8-40 Mq (Heger et al. 2003). 
The stars with mass below or above this mass range become 
AGB stars or 'coUapsers', respectively. 

We neglect Type la SNe for simpUcity. This population 
of SNe is the major source of iron element (Iwamoto et 
al. 1999) and may be the source of iron dust (Calura et 
al. 2008). However, in respect of the total Stardust mass 
budget, the contribution relative to SNe II is always less than 
1-10% (Zhukovska et al. 2008; Pipino et al. 201 1). Since we 
are dealing with metal and dust as each a single component, 
we safely neglect the contribution of SNe la. 

The remnant mass, w(m, Z), is taken from model calcu- 
lations of AGB stars (Karakas 2010) and SNe (Nomoto et 
al. 2006). Figure 1 shows the remnant mass fraction relative 
to the initial stellar mass, w/m. This depends on metallicity 
Z because Z in the stellar atmosphere determines radiation 
pressure through opacity and the strength of the stellar wind 
in the course of the stellar evolution, and affects the rem- 
nant mass. However, as shown in Figure 1, the dependence 
is weak, so we neglect it. We obtain the following fitting 
formula: 



iu{m) 



1 

0.13 



m 



8M(7 







1 •? ( ™ 



-0.5 



-0.7 



(m > 40 Mq) 

(8 Mq < m < 40 Mq) 

(m < 8 Mq) 



(19) 

which is shown by the solid line in Figure 1. This fitting 
formula agrees with the values in Table 1 of Morgan & 
Edmunds (2003) within a < 20% difference, except for m = 
9 Mq case in which our estimate is a factor of 2 lower than 
that of Morgan & Edmunds (2003). 

For the metal yield, mz, we adopt the data taken 
from model calculations of AGBs (Karakas 2010) and SNe 
(Nomoto et al. 2006). Figure 2 shows mz relative to the ini- 
tial stellar mass m as a function of m. While the expected 
mz depends on mass m and metallicity Z by a complex way. 
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Fig. 1. Remnant mass fraction, w/m, as a function of the initial stellar 
mass, m. The data of AGB stars (1-8 Mq) are taken from Karakas 
(2010) and those of SNe (8^0 Mq) are taken from Nomoto et al. (2006). 
The different symbols indicate different metallicity Z. Z = (plus), 
Z = 0.0001 (cross), Z = 0.001 (diamond), Z = 0.004 (triangle), 
Z = 0.008 (circle), and Z = 0.02 (square). The solid line is a fitting 
function given by equation (19). 




1 10 100 

Initial stellar mass [Mq\ 



Fig. 2. Ejected metal mass fraction, mz/m, as a function of the initial 
stellar mass, m. The data of AGB stars (1-8 Mq) are taken from Karakas 
(2010) and those of SNe (8^0 Mq) are taken from Nomoto et al. (2006). 
The meaning of the symbols are the same as in Fig. 1. The solid line is 
a fitting function given by equation (20). The dotted lines are the cases a 
factor of two higher or lower than the solid line. 

we approximate the data with a simple power-law of only m 

as 

(0 (m > 40 Mq) 

/^(s^)' (8M0<m<4OM0) 




(m < 8Mo) 



(20) 

When the normalization fz — 0.02, equation (20) is the 
solid line in Figure 2. As shown in the figure, the uncertainty 
of equation (20) is a factor of ^ 2. This fitting agrees with 



Fig. 3. Ejected dust mass fraction, /m, as a function of the initial stellar 
mass, m. The data of AGB stars (1-8 Mq) are taken from Zhukovska 
et al. (2008) and those of SNe (8^0 Mq) are taken from Nozawa et 
al. (2007). For the AGB data, the different symbols indicate different 
metaUicity Z: Z = 0.001 (diamond), Z = 0.004 (triangle), Z = 0.008 
(circle), and Z = 0.02 (square). For the SNe data, the different symbols 
indicate different ambient hydrogen density nu = 0.1 cm""^ (circle), 
njj = 1 cm~^ (triangle), and rtu = 10 cm~^ (square). The open and 
filled symbols correspond to 'mixed' or 'unmixed' cases of Nozawa et 
al. (2007), respectively. The dotted, solid, and dashed lines are the cases 
of ^ = 1, 0.1, and 0.001, respectively, in equation (21). 

the values in Table 1 of Morgan & Edmunds (2003) within 
a factor of 2 difference in the SN regime. However, in the 
AGB regime, the difference is as large as the model results 
by Karakas (2010). The effect of this large uncertainty of the 
yield is discussed in §6.1. 

The dust yield, rrid, calculated by Zhukovska et al. (2008) 
and Ferrarotti & Gail (2006) for AGBs and Nozawa et 
al. (2007) for SNe are shown in Figure 3. These yields are 
theoretical ones and do not seem to be compared with obser- 
vations very much yet. As found in Figure 3, depends 
on mass m and metallicity Z by a complex way as the metal 
yield mz does. Moreover, the dust production by SNe is fur- 
ther complex because the reverse shock moving in the ejecta 
of a SN may destroy the dust produced in the ejecta (Bianchi 
& Schneider 2007, Nozawa et al. 2007, Nath et al. 2008, 
Silvia et al. 2010). This self-destruction depends on the ma- 
terial strength against the destructiord and the ambient gas 
density which determines the strength of the reverse shock. 
According to Nozawa et al. (2007), we plot three cases of 
the ambient density and 'mixed' and 'unmixed' dustproduc- 
tionflin Figure 3. We adopt a simple formula for as 

md{m) = £,mz{m) , (21) 

where ^ is a scaling factor and means an efficiency of con- 



^The micro-process of the destruction considered in Nozawa et al. (2007) is 
sputtering by hot gas. 

^The 'mixed' and 'unmixed' mean the elemental mixing in the SN ejecta 
(Nozawa et al. 2003). In the 'mixed' case, there is no layer where C is more 
abundant than O, then, only silicate, troilite, and corundum grains can be 
formed. On the other hand, the 'unmixed' case has a C-rich layer and Fe 
layer and can form carbon and iron grains as well as silicate. 



6 



A. K. INOUE: THE ORIGIN OF DUST IN GALAXIES 




Fig. 4. Stai' formation history (time evolution of star formation rate) 
normalized by the average rate. The histogram is the observed history 
at the solar neighborhood reported by Rocha-Pinto et al. (2000a). The 
filled circles with eiTor-bars are the average of the histogram over 4 Gyr 
interval and corrected by the average rate only for the age less than 12 
Gyr which is the present age of the Milky Way assumed in this paper. 
The vertical error-bars are the standard error of the mean. The four lines 
correspond to the time evolutions with four different sets of time-scales 
of star formation and infall (rsp/Gyr, Tin/Gyr) as indicated in the panel. 

densation of metal elements. In Figure 3, we show three 
cases of ^ = 1 (all metal condenses into dust; an extreme but 
unrealistic case), 0.1 (fiducial case), and 0.01 (a lower effi- 
ciency case). The reader may be anxious about a large un- 
certainty of this approximation. However, the dust mass in 
galaxies does not depend on nid after the accretion growth 
becomes active. This is because the growth of dust is the 
dominant process of dust production after the activation as 
shown later in §6.1. 

4. Milky Way analog 

Let us calibrate parameters in the chemical and dust evolu- 
tion model of galaxies so as to reproduce the properties at the 
solar neighborhood in the Milky Way. There are two param- 
eters in the chemical evolution part: the time-scales of star 
formation, rgp, and infall. Tin. There are additional two pa- 
rameters in the dust content evolution: the time-scale of the 
ISM accretion growth, Tac,Oi and the efficiency of the dust 
destruction, etosn- In addition, there are two parameters as 
uncertainties of metal and dust yields, fz and ^. Table 1 is a 
summary of these parameters and values. 

Note that we do not apply any statistical method to justify 
the goodness of the reproduction of the observational con- 
straints throughout this paper because our aim is not to find 
the best fit solution for the constraints but to demonstrate the 
dust content evolution in galaxies qualitatively. This is partly 
due to the weakness of the observational constraints and due 
to large uncertainties of the dust physics itself. 
4.1 Chemical evolution at the solar neighborhood 

We here determine the time-scales of star formation and 
infall in the chemical evolution part. First, we constrain 
these time-scales by using the star formation history at the 
solar neighborhood reported by Rocha-Pinto et al. (2000a). 



Fig. 5. Age-metallicity relation of stars. The diamonds with error-bars are 
the data of stars at the solar neighborhood reported by Rocha-Pinto et 
al. (2000b). The solar mark (0) indicates the position of the Sun on this 
plot. The four lines correspond to the model relations same as Fig. 4. The 
present age of the Milky Way is assumed to be 1 2 Gyr. 



Such a method was adopted by Takeuchi & Hirashita (2000). 
Rocha-Pinto et al. (2000a) derived the star formation history 
from the age distribution of 552 late-type dwarf stars at the 
solar neighborhood. The histogram in Figure 4 is their result 
and shows very stochastic nature of the history. However, our 
model can treat only a smooth history. Thus, we smoothed 
the stochastic history by averaging with 4 Gyr interval. The 
filled circles are the result. The vertical error-bars indicate 
the standard error of the mean. The average history is re- 
normalized by the average star formation rate for the stellar 
age less than 12 Gyr which is the assumed age of the Milky 
Way in this paper, although this choice of the age is arbitrary. 
We have tried four cases of tsf in this paper: 1, 2, 3, and 
5 Gyr which are the observed range of the time-scale (or 
gas consumption time-scale) for disk galaxies like the Milky 
Way (e.g., Larson, Tinsley, & Caldwell 1980). For each 
Tsf, we have found Tin with which we can reproduce the 
smoothed history as shown in Figure 4. 

Next, we adopt the observed relation between the stellar 
age and metallicity, so-called the age-metallicity relation, 
reported by Rocha-Pinto et al. (2000b) to further constrain 
(tsf, Tin). Rocha-Pinto et al. (2000b) derived the relation 
from the same 552 stars as Rocha-Pinto et al. (2000a). Their 
result is shown in Figure 5 by diamonds with error-bars. 
After comparing with our four model lines, we have found 
that (TgF/Gyr, Tin/Gyr) = (5,10) case seems the best 
match with the observed relation but (tsf / Gyr, Tin/Gyr) — 
(3, 15) case is also acceptable. 

Finally, we adopt another constraint: the current stellar 
mass relative to the ISM mass. Naab & Ostriker (2006) com- 
piled observational constraints for the solar neighborhood. 
From the compilation, we adopt the ratio of the stellar mass 
to the ISM mass at the present epoch of 3.1 ± 1.4. Note that 
the stellar mass includes the remnant mass (i.e. white dwarf, 
neutron stars, and black-holes). Figure 6 shows the com- 
parison of the ratio with our four star formation histories. 
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Table 1 . Parameters and values for the solar neighborhood. 



Parameter 


Fiducial value 


Considered values 


(rsp/Gyr, rin/Gyr) 


(3, 15) 


(1,50), (2, 20), (3, 15), and (5, 10) 


(Tac,o/Myr, emsN/lO^M©) 


(3, 1) 


(1.5, 1), (1.5, 2), (3, 0.5), (3, 1), (3, 2), (6, 0.5), and (6, 1) 


fz 


0.02 


0.01, 0.02, and 0.04 




0.1 


0.01, 0.1, and 1 




Fig. 6. Time evolution of the stellar mass relative to the ISM mass. The 
filled circle with error-bars is an estimate at the solar neighborhood taken 
from Naab & Ostriker (2006) and references therein. Note that the stellar 
mass includes mass of remnants. The present age of the Milky Way 
should be less than the age of the Universe (13.7 Gyr). The four lines 
con'espond to the model evolutions same as Fig. 4. 



We have found that the two sets of (rsp/Gyr, Tin/Gyr) = 
(2, 20) and (3,15) are consistent with the data. 

From these three comparisons, we finally adopt the case 
of (rsp/Gyr, Tin/Gyr) = (3, 15) as the fiducial set for the 
Milky Way (or more precisely, for the solar neighborhood) 
in this paper 

4.2 Dust content evolution at the solar neighborhood 

Here we examine the dust content evolution. First, we 
show the significant effect of the dust destruction and the 
ISM growth. Figure 7 shows the time evolution of metallicity 
and dust-to-gas mass ratio for the fiducial set of tsf and 
Tin obtained in the previous subsection. The model curves 
of the dust-to-gas ratio (dotted, dot-dashed, and solid lines) 
are compared with the filled circle with error-bars which is 
an observational estimate at the solar neighborhood. This is 
obtained from metallicity Z k, Zq (van den Bergh 2000; see 
also Rocha-Pinto et al. 2000b) and dust-to-metal mass ratio 
(5 w 0.5 (Kimura et al. 2003; see below) and the uncertainty 
is the quadrature of uncertainties of 309^ in Z and 20% in 
S. 

If there is neither destruction nor accretion growth of dust, 
the dust-to-gas ratio evolution is just the metallicity evolution 
multiplied by the condensation efficiency of Stardust, ^, as 



^The difference between Zqs by Anders & Grevesse (1989) and Asplund 
et al. (2009) accounts for the uncertainty. 



Fig. 7. Time evolution of metallicity (metal mass fraction in the ISM; 
dashed fine) and dust-to-gas mass ratio (other lines) for the fiducial case 
(Table 1). The dotted fine corresponds to the case only with 'Stardust' 
production and destruction by star formation (i.e. astration). The dot- 
dashed line correspond to the case with the dust destruction by SNe but 
without the dust growth in the ISM. The solid fine is the case with all the 
processes. The filled circle with eiTor-bars is an estimate of the dust-to- 
gas mass ratio at the solar neighborhood (see text). The present age of 
the Milky Way should be less than the age of the Universe (13.7 Gyr). 



shown by the dotted line. We have assumed ^ = 0.1 for the 
line. Once the SN destruction of dust is turned on with a 
standard efficiency as emsN = 1 x 10'^ Mq (McKee 1989; 
Nozawa et al. 2006), it reduces the dust amount by a factor of 
ten as shown by the dot-dashed line. This confirms that the 
dust destruction is very efficient and the Stardust injection is 
too small to compensate the destruction (e.g., Draine 1990; 
Tielens 1998). Then we need the accretion growth in the 
ISM to reproduce the dust-to-gas ratio 10"^ in the present 
Milky Way. If we assume the time-scale of Tac.o = 3 x 10^ 
yr, the dust-to-gas ratio evolution becomes the solid line and 
it reaches ~ 10~^ which is almost two orders of magnitude 
larger than the case without the growth after several Gyr 

Figure 8 shows the time evolution of tsn in equation (13) 
and Tac in equation (17). The SN destruction time-scale tsn 
is almost constant promptly after the first a few hundreds 
Myr On the other hand, the accretion growth time-scale Tac 
decreases gradually in the first a few Gyr. This is because 
Tac has a metallicity dependence as shown in equation (17) 
and it decreases as the metallicity increases. At the time 
around 4 Gyr, Tac becomes shorter than tsn, and then, the 
accretion growth becomes significant and the dust amount 
increases rapidly. As the accretion growth proceeds, the 
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Fig. 8. Time evolution of time-scales of tiie dust destruction by SNe (daslied 
line) and of tlie dust growth in the ISM (solid line) for the fiducial case 
(Table 1). 




Fig. 9. Time evolution of dust-to-metal mass ratio, 5, for various sets of 
the parameters of dust growth and destruction (rac.o, fiSN)- The solid 
line is the fiducial case with (3 X 10'' yr, 1 X lO^ Mq). The two dotted 
lines con'espond to two different sets of (1.5 X lO'' yr, 1 X 10'^ Mq) 
for the upper and (1.5 X 10*' yr, 2 X 10^ Mq) for the lower. The two 
dot-dashed lines correspond to the sets of (3 X 10® yr, 5 X lO'^ Mq) 
for the upper and (3 X 10'^ yr, 2 X 10^ Mq) for the lower. The two 
dashed lines correspond to the sets of (6 X 10® yr, 5 X 10^ Mq) for the 
upper and (6 X 10® yr, 1 X 10^ Mq) for the lower. The filled circle with 
error-bars is an estimate in the Local Interstellar Cloud surrounding the 
Sun by Kimura et al. (2003). The present age of the Milky Way should 
be less than the age of the Universe (13.7 Gyr). 



metal abundance in the gas phase decreases, i.e., the dust- 
to-metal ratio 5 increases, then, Tac becomes almost constant 
and balances with tsn. We will discuss this point in §5 more 
in detail. 

Figure 9 shows the time evolution of the dust-to-metal ra- 
tio, 6. The solid line is the fiducial case which is shown in 
Figures 7 and 8. This can be compared with the observed 
ratio in the Local Interstellar Cloud reported by Kimura et 
al. (2003): (5 = 0.5 ± 0.1. As shown in Figure 9, the fidu- 



cial set of (Tac,o, e"isN) = (3 Myr, 1 x 10'^ Mq) excellently 
agrees with the observed data. On the other hand, other sets 
can also reproduce the data. For example, (racOi e^sn) = 
(1.5 Myr, 2 x IO^A/q) or (6 Myr, 5 x IO^A/q). Interest- 
ingly, the 5 evolutions become very similar if the product of 
Tac and emsN is the same. We will also discuss this point 
in §5. 

5. Determining dust-to-metal ratio 

In this section, we demonstrate the mechanism for deter- 
mining the dust-to-metal mass ratio, 5, in galaxies. Starting 
from equations (3) and (4), we can obtain the time evolution- 
ary equation of 5 = M^jMz- Here, let us adopt the instan- 
taneous recycling approximation (IRA) in which we neglect 
the finite stellar life-time and assume that stars with a mass 
larger than a certain threshold (the turn-off mass rrit) die in- 
stantly when they are formed. This approximation allows 
us to manage the equations analytically and is good enough 
to see phenomena with a time-scale longer than a Gyr (see 
Tinsley 1980 for more details). In the IRA, we can approxi- 
mate the metal and dust yields in equations (10) and (11) as 
Yz ~yzS and = iYz « ^yzS, where the effective 
metal yield 



yz 



mz(rn)(f>(m)d'm = 0.024 



fz 



0.02 



(22) 



where we have assumed mt = I MQ. This value is not 
sensitive to We obtain yz = 0.021 (/2/O.O2) if = 5 
Mq. Remembering the star formation rate 5* ~ Mism/tsf 
as in equation (5), then, we obtain 



ldS_ 
Sdt 



yz 

TspZ 



-I 



1 

TSN 



-■ (23) 



In the IRA, the SN destruction time-scale tsn in equation 
(13) can be reduced to 



TSF 



emsN^sN 



(24) 



where the effective number of SN per unit stellar mass is 

•4OM0 

'SMq 



/■4OM0 

/ (j){m)dm = 0.010 Mq-^ . (25) 

J 8 Mm 



Note that emgN^^sN is a non-dimensional value. The accre- 
tion growth time-scale Tac is given in equation (17). Then, 
equation (23) is reduced to 



ld5_ 
Sdi 



a + emsN^sN , Z{1 - 6) 



where 



TSF 

yz 
z 



Tac.O 



1 - 



(26) 



(27) 



In the IRA, the metallicity Z = A/^/A/jsm can be ob- 
tained analytically (for example, see Dwek et al. 2007). 
Then, we have found that Z — > yz for < — !■ 00 when 
Tin > Tsp. The condensation efficiency ^ is uncertain but 
it is of the order of 0. 1 (see Figure 3). When S is of the order 
of 0. 1-1 as shown in Figure 9, the ratio ^/(5 is of the order of 
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1 or smaller. Therefore, a is also of the order of 1 or smaller. 
On the other hand, nsN ^ 10^^ M^ '^ and emsN ^ 10^ M( 



'0. 



then, we obtain emsN?^SN S> oi. Therefore, equation (26) is 
further reduced to 



ld5_ 
Sdi 



+ 6(1-5) 



(28) 



where a = etosn't^sn/tsf and b = Z/t^co- If we assume Z 
to be constant (i.e. b is constant), equation (28) can be solved 
analytically. The solution is 



SooSq cxp {b — a)t 
{5oo - 6q) + 6o exp {b - a)t 



(29) 



where 5q and 5oo are the values for < = and t — > cx), 
respectively. The asymptotic value 5oo for t — !■ cx) is realized 
only when b > a, and is given by 



l-S^ = - = 



T■ac,oe'7^sN?^SN 



(30) 



b tsfZ 
This is the equilibrium value for equation (28) and we find 



1 



0.5 



3 Myr J 



emsN 



{iO^MqJ ^10 
3Gyr 

TSF 



Me 
0.02 



(31) 



which excellently agrees with the results in Figure 9. 

We can fully understand the S evolution by using equation 
(28). At the beginning, the accretion term 6^0 because 
Z ~ 0. Then, only the destruction term a is effective. As a 
result, S decreases with the time-scale of 1/a — tsn- As Z 
increases, the accretion term b increases and finally exceeds 
a. Then, 6 increases toward Soo with the evolution time- 
scale of 1/(6 — a). This decreases as Z increases and 6 — a 
increases. Therefore, the driving force of the S evolution is 
Z. If we call Z at b = a as the critical metallicity, Zc, we 
find 



;,oemsN'T-SN 



= 0.01 




When Z > Zc, the accretion growth becomes effective and S 
approaches the final value Soo- A similar critical metallicity 
has been derived by Asano et al. (201 1) with a different way. 

Equation (30) shows that the final value of S is determined 
by the equilibrium between the SN destruction and the ac- 
cretion growth in the ISM. The time-scale to reach the equi- 
librium is 1/(6 — a). This is relatively short in the fiducial 
case. For example, it is 0.3 Gyr when Z = 0.02. This 
means that the 6 evolution proceeds with keeping the equilib- 
rium between the SN destruction and the accretion growth, 
or equivalently, 5 — Soo after Z exceeds Zc- This behavior 
is also found by the comparison of the two time-scales, tsn 
and Tac, in Figure 8; once Tac becomes shorter than tsn at 
about 4 Gyr at which Z exceeds Zc, Tac turns around and ap- 
proaches TSN again. This is realized by the reduction of the 
term (1 — 5) in Tac (see eq. [17]) when S increases from ^ 
to Soo - Such a kind of self -regulation process determines the 
dust-to-metal ratio S. 



. ' ' ' I ' ' ' I ' ' ' I ' ' 

; (a) (Tsr,T,„) = (3,,l-St' 




4 6 8 10 
Time [Gyr] 



Fig. 10. Time evolution of (a) metallicity, (b) dust-to-metal mass ratio, and 
(c) time-scales of SN destruction and ISM growth for three different stel- 
lar metal yields: fz = 0.01 (dot-dashed), 0.02 (solid), and 0.04 (dotted). 
The horizontal short-dashed line in the panel (a) shows the critical metal- 
Ucity of equation (32). The vertical short-dashed lines indicate the timing 
at which the metalUcity exceeds the critical one for the three metal yields. 
The long-dashed line in the panel (c) is the SN destruction time-scale, but 
the other three lines are the ISM growth time-scales for the three metal 
yields. 



6. Discussion 

6.1 Effect of uncertainties of yields 

Here we examine the effect of uncertainties of the normal- 
ization of metal and dust yields. As we saw in Figure 2, our 
simple recipe for the metal yield may contain a factor of 2 (or 
more) uncertainty. The parameter fz accounts for this uncer- 
tainty. In Figure 10, we show the effect of fz- As found from 
the panel (a), the metallicity evolution is scaled almost lin- 
eally by fz as expected and the timing at which Z exceeds 
Zc = 0-5Zq given by equation (32) for the fiducial set of 
the accretion and destruction efficiencies becomes faster as 
fz is larger From the panels (b) and (c), we find that for 
each case of fz, S increases and Tac becomes shorter than 
TgN soon after the timing for Z > Zc- Therefore, the timing 
for Tac < TSN, in other words, the timing for the accretion 
growth activation is well traced by Zc in equation (32) and 
this is not affected by uncertainty of fz- On the other hand, 
the timing for the activation becomes faster for larger fz- 
The metallicity dependence on the final value of 5 is explicit 
as found in equation (31). 

Figure 1 1 shows the effect of the dust yield. As seen in 
Figure 3, our recipe for the Stardust yield has a factor of 10 
or larger uncertainty because of a large uncertainty in the 
adopted model calculations. In Figure 1 1, we show the cases 
with a factor of 10 larger or smaller yield than the fiducial 
one. Other parameters are the same as the fiducial set, so 
that we have the same evolutions of the metallicity and the 
time-scales of the SN destruction and accretion growth as 
shown by the solid lines in Figure 10. Before the growth 
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Fig. 11. Same as Fig. 7 but for different condensation efficiencies in the 
stellar ejecta as indicated in the panel. 



activation at around 4 Gyr, the dust amounts sho'w a large 
difference, ho'wever, they converge nearly the same amount 
after the activation. This is because the final value of S 
given in equation (31) does not depend on the dust yield. 
Therefore, we conclude that the dust content in galaxies is 
independent of the Stardust yield after the grain growth in the 
ISM becomes active, or equivalently, the metallicity exceeds 
the critical one. 

6.2 What kind of dust is formed by the ISM growth? 

We have shown that the main production channel of dust 
is the accretion growth in the ISM of the present-day Milky 
Way. This conclusion had been obtained also in the litera- 
ture. For example, Zhukovska et al. (2008) argued that the 
mass fraction of stardusts in total dust is only 0.1-1% based 
on a more sophisticated chemical evolution model than this 
paper (see their Fig. 15); more than 99% of dust is originated 
from the accretion growth in the ISM. It is also well known 
that some interplanetary dust particles show a highly en- 
hanced abundance of deuterium and ^^N relative to the solar 
composition, which is a signature of the molecular cloud ori- 
gin because such isotopic fractionations are expected in low 
temperature environment (e.g.. Messenger 2000). Therefore, 
dust produced by the ISM accretion exists. Then, we have 
a very important question; what kinds of dust species are 
formed by the accretion growth in the ISM? 

In molecular clouds, many kinds of ices such as H2O, CO, 
CO2, CH3OH have been detected (e.g., Gibb et al. 2000). 
These ices are condensed onto pre-existent grains. In these 
ices, some chemical reactions and ultraviolet photolysis (and 
cosmic rays) process the material and may make it refrac- 
tory. As a result, so-called 'core-mantle grains' coated by 
refractory organics would be formed (e.g., Li & Greenberg 
1997). Indeed, such a grain has been found in cometary dust: 
olivine particles produced by a Type II SN coated by organic 
matter which seems to be formed in a cold molecular cloud 
(Messenger et al. 2005). Therefore, the ISM dust probably 
has core-mantle or layered structures. Moreover, the com- 
position can be heterogeneous: for example, graphite coated 



by silicate, silicate coated by graphite, silicate coated by iron, 
etc. The formation of such grains does not seem to be stud- 
ied well. Much more experimental and theoretical works are 
highly encouraged. 

If we can find signatures of the dust accretion growth in 
the ISM of galaxies by astronomical observations (i.e. very 
distant remote-sensing), it proves the growth ubiquitous. A 
possible evidence already obtained is a huge mass of dust in 
galaxies which requires the accretion growth as discussed in 
this paper It is worth studying how to distinguish Stardust 
grains (or grain cores) and ISM dust (or mantle) by observa- 
tions, e.g., spectropolarimetry, in future. 
6.3 Dust amount in the proto-solar nebula 

We have shown that the dust amount is very small before 
the ISM growth becomes active. For example, the dust-to- 
gas mass ratio is of the order of lO^"* at the first a few Gyr 
from the formation of the Milky Way (or the onset of the 
major star formation at the solar neighborhood). If the dust- 
to-gas ratio in the proto-solar nebula was 10^"*, the planet 
formation might be difficult. Fortunately, the activation of 
the ISM growth is expected to be about 8 Gyr ago at the solar 
neighborhood. Thus, it is well before the solar system forma- 
tion. Indeed, we expect the dust-to-gas ratio of several times 
lO^'^ at 4-5 Gyr ago (see Figure 7). Moreover, the dust-to- 
gas ratio may be much enhanced in the proto-solar nebula 
relative to the average ISM. This is because the accretion 
growth is more efficient for higher density and the density in 
the proto-solar nebula is several orders of magnitude higher 
than that in molecular clouds. Therefore, even if the solar 
system formation is before the activation of the ISM growth 
globally, the dust growth may be active locally in the proto- 
solar nebula. In this case, the planet formation is always pos- 
sible if there is enough metal to accrete onto the pre-existent 
seed grains, even before the global growth activation. This is 
an interesting issue to relate to the Galactic Habitable Zone 
where complex life can be formed (Lineweaver et al. 2004). 
We will investigate it more in future. 
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